Differential haplotype expression in class I MHC genes during SARS-CoV-2 infection of human lung cell lines (preprint)

Cell entry of SARS-CoV-2 causes genome-wide disruption of the transcriptional profiles of genes and biological pathways involved in the pathogenesis of COVID-19. Expression allelic imbalance is characterized by a deviation from the Mendelian expected 1:1 expression ratio and is an important source of allele-specific heterogeneity. Expression allelic imbalance can be measured by allele-specific expression analysis (ASE) across heterozygous informative expressed single nucleotide variants (eSNVs). ASE reflects many regulatory biological phenomena that can be assessed by combining genome and transcriptome information. ASE contributes to the interindividual variability associated with disease. We aim to estimate the transcriptome-wide impact of SARS-CoV-2 infection by analyzing eSNVs. We compared ASE profiles in the human lung cell lines Calu-3, A459, and H522 before and after infection with SARS-CoV-2 using RNA-Seq experiments. We identified 34 differential ASE (DASE) sites in 13 genes (HLA-A, HLA-B, HLA-C, BRD2, EHD2, GFM2, GSPT1, HAVCR1, MAT2A, NQO2, SUPT6H, TNFRSF11A, UMPS), all of which are enriched in protein binding functions and play a role in COVID-19. Most DASE sites were assigned to the MHC class I locus and were predominantly upregulated upon infection. DASE sites in the MHC class I locus also occur in iPSC-derived airway epithelium basal cells infected with SARS-CoV-2. Using an RNA-Seq haplotype reconstruction approach, we found DASE sites and adjacent eSNVs in phase (i.e., predicted on the same DNA strand), demonstrating differential haplotype expression upon infection. We found a bias towards the expression of the HLA alleles with a higher binding affinity to SARS-CoV-2 epitopes. Independent of gene expression compensation, SARS-CoV-2 infection of human lung cell lines induces transcriptional allelic switching at the MHC loci. This suggests a response mechanism to SARS-CoV-2 infection that swaps HLA alleles with poor epitope binding affinity, an expectation supported by publicly available proteome data.

Genomic perspectives on the emerging SARS-CoV-2 Omicron variant (preprint)

A new variant of concern for SARS-CoV-2, Omicron (B.1.1.529), was designated by the World Health Organization on November 26, 2021. This study analyzed the viral genome sequencing data of 108 samples collected from patients infected with Omicron. First, we found that the enrichment efficiency of viral nucleic acids was reduced due to mutations in the region where the primers anneal to. Second, the Omicron variant possesses an excessive number of mutations compared to other variants circulating at the same time (62 vs. 45), especially in the Spike gene. Mutations in the Spike gene confer alterations in 32 amino acid residues, which was more than those observed in other SARS-CoV-2 variants. Moreover, a large number of nonsynonymous mutations occur in the codons for the amino acid residues located on the surface of the Spike protein, which could potentially affect the replication, infectivity, and antigenicity of SARS-CoV-2. Third, there are 53 mutations between the Omicron variant and its closest sequences available in public databases. Many of those mutations were rarely observed in the public database and had a low mutation rate. In addition, the linkage disequilibrium between these mutations was low, with a limited number of mutations (6) concurrently observed in the same genome, suggesting that the Omicron variant would be in a different evolutionary branch from the currently prevalent variants. To improve our ability to detect and track the source of new variants rapidly, it is imperative to further strengthen genomic surveillance and data sharing globally in a timely manner.

Pervasive transmission of E484K and emergence of VUI-NP13L with evidence of SARS-CoV-2 co-infection events by two different lineages in Rio Grande do Sul, Brazil (preprint)

Emergence of novel SARS-CoV-2 lineages are under the spotlight of the media, scientific community and governments. Recent reports of novel variants in the United Kingdom, South Africa and Brazil (B.1.1.28-E484K) have raised intense interest because of a possible higher transmission rate or resistance to the novel vaccines. Nevertheless, the spread of B.1.1.28 (E484K) and other variants in Brazil is still unknown. In this work, we investigated the population structure and genomic complexity of SARS-CoV-2 in Rio Grande do Sul, the southernmost state in Brazil. Most samples sequenced belonged to the B.1.1.28 (E484K) lineage, demonstrating its widespread dispersion. We were the first to identify two independent events of co-infection caused by the occurrence of B.1.1.28 (E484K) with either B.1.1.248 or B.1.91 lineages. Also, clustering analysis revealed the occurrence of a novel cluster of samples circulating in the state (named VUI-NP13L) characterized by 12 lineage-defining mutations. In light of the evidence for E484K dispersion, co-infection and emergence of VUI-NP13L in Rio Grande do Sul, we reaffirm the importance of establishing strict and effective social distancing measures to counter the spread of potentially more hazardous SARS-CoV-2 strains. HighlightsO_LIThe novel variant B.1.1.28 (E484K) previously described in Rio de Janeiro is currently spread across the southernmost state of Brazil; C_LIO_LIThe novel variant VUI-NP13L was also identified by causing a local outbreak in Rio Grande do Sul; C_LIO_LIB.1.1.28 (E484K) is able to establish successful coinfection events co-occurring simultaneously with different lineages of SARS-CoV-2. C_LI

Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil (preprint)

In this study, we report the sequencing of 180 new viral genomes obtained from different municipalities of the state of Rio de Janeiro from April to December 2020. We identified a novel lineage of SARS-CoV-2, originated from B.1.1.28, distinguished by five single-nucleotide variants (SNVs): C100U, C28253U, G28628U, G28975U, and C29754U. The SNV G23012A (E484K), in the receptor-binding domain of Spike protein, was widely spread across the samples. This mutation was previously associated with escape from neutralizing antibodies against SARS-CoV-2. This novel lineage emerged in late July being first detected by us in late October and still mainly restricted to the capital of the state. However, as observed for other strains it can be rapidly spread in the state. The significant increase in the frequency of this lineage raises concerns about public health management and continuous need for genomic surveillance during the second wave of infections.